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Abstract. Compared to planets around Sun-like stars, relatively little is 
known about the occurrence rate and orbital properties of planets around stars 
more massive than 1.3 Mq. The apparent deficit of planets around massive stars 
is due to a strong selection bias against early-type dwarfs in Doppler-based 
planet searches. One method to circumvent the difficulties inherent to massive 
main-sequence stars is to instead observe them after they have evolved onto 
the subgiant branch. We show how the cooler atmospheres and slower rotation 
velocities of subgiants make them ideal proxies for F- and A-type stars. We 
present the early results from our planet search that reveal a paucity of planets 
orbiting within 1 AU of stars more massive than 1.5 Mq, and evidence of a rising 
trend in giant planet occurrence with stellar mass. 



Introduction 



A planet-bearing star can be thought of as a very bright, extremely dense rem- 
nant of a protoplanetary disk. After all, a star inherits its defining characteristics- 
its mass and chemical composition — from the same disk material that forms its 
planets. The physical characteristics of planet host stars therefore provide a 
crucial link between the planets we detect today and the circumstellar environ- 
ments from which they formed long ago. Studying the relationships between 
the observed occurrence rate and orbital properties of planets as a function of 
stellar characteristics informs theories of planet formation, and also helps guide 
the target selection of future planet searches. 

A wealth of recent work has demo nstrated that planet occurrence is strongly 
correlated with chemical composition ( Gonzale^ll997l : ISantos et al.ll2004l ): metal- 
rich stars are 3 times more likely to host planetary companions compared to 
stars with solar abundances ( Fischer &: Valenti 20051 ). This finding can be un- 
derstood in the context of the core accretion model. Increasing the metallicity 
of a star/disk system increases the surface density of solid material at the disk 
midpl ane, which in turn leads to an enh anced growth rate for protoplanetary 
cores (|lda Linll2"oo3 : iKornet et al.ll2005l ^. 

Another factor that enhances the surface density of material in the disk 
midplane is its total mass. If the mass of circumstellar disks scales with the 
mass of the central star, then there sh ould be an observed correlation be- 
tween planet occurrence a nd stellar mass (iLaughlin et al.ll2004l : llda Lin|[2005l : 
Kennedy &: Kenvon 20071 ). In principle, testing this hypothesis is fairly simple: 
one need only measure the fraction of stars with planets over a wide range of 
stellar masses. However, in practice such a study is not so straight forward given 
the limited range of stellar masses encompassed by most planet searches. 
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The difficulty can be seen in Figure [U which shows the distribution of 
stehar r nasses for the target st ars in Cahfornia and Carnegie Planet Search 
(CCPS; IValenti Sz Fischer |2005| ). which is representative of most Doppler-based 
planet searches. Most of the stars in Figure [T] have masses between 0.7 Mq and 
1.3 Mq. In a decidedly non-Copernican twist of nature, it turns out that stars 
like our Sun are ideal planet search targets. Solar-mass G and K dwarfs are 
slow rotators, have stable atmospheres, and are relatively bright. The fall-off 
toward lower stellar masses is simply because late K- and M-type dwarfs are 
faint, making the a cquisition of h igh-resolution spectra d ifficult without l arge 
telescope apertures (iButler et al.l [2006: B onfils et al. 2005.: lEndl et al.l 120031 '). 

The sharp drop at higher stellar masses is due to a separate observational 
bias. Stars with spectral types earlier than F8 tend to have rotationally broad- 
ened spectral features {Vsini > 50 km s~^; Galland et al. 2005), have fewer 
spectral lines due to high surface temperatures, and display a large amount of 
atmospheric "jitter." Stellar jitter is excess velocity scatter due to surface inho- 
mogeneities and pulsation, whi ch can approach 50-100 m for mid-F stars 
(|Saar et al.l Il998l : IWrightl l2005l ). These features conspire to limit the attain- 
able radial velocity precision for early-type dwarfs to > 50 m s~^, rendering 
exceedingly difficult the detection of all but the most massive and short-period 
planets. 

One method to circumvent the observational limitations inherent to high- 
mass dwarfs is to observe these stars after they have evolved away from the 
main-sequence. After stars have expended their core hydrogen fuel sources their 
radii expand, and their atmospheres cool leading to an increase in the number 
of metal lines in the star's spectrum. Stars crossing the subgiant branch also 
shed a large amount of angular momentur n through the coupling of stellar winds 
to rotationally generated inagne tic fields ( Grav &: NagaijlT985l : Schrijver &: Fold 
19931 : do Nascimento et al. 2000l ). The cooler atmospheres and slower rotational 
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velocities of evolved stars lead to an increased number of narrow absorption lines 
in their spectra, making them much better suited for precision Doppler surveys. 



2. A Doppler Survey of Subgiants 



There are a number of planet searches targeting evolved, 
stars. To date, niost surveys have focused on K-gian ts 



Liu et al 



inter mediate-mass 
( Frink et al l l2002l : 



Hatzes et alll2005l : iLovis fc Mavoill2007l : iNiedzielski et al.ll2007l) and "clump gi- 
ants," or asym ptotic giant branch stars ( Sato et al. 20031 : Setiawan et al.l [20031 : 



20071 1. These programs have detected a total of 9 substellar compan- 
ions orbiting intermediate-mass giants, demonstrating thatplanets do form and 
can be detected around stars more massive than ~ 1.5 Mgu- 

Over the past 3 years we have been conducting a Doppler survey of evolved 
stars at Lick and Keck Observatories. However, instead of targeting clump giants 
and K giants, we have focused on stars occupying the region of the H-R diagram 
between the main-sequence and red giant branch , known as the subgiant branch. 
Our sample is described bv I Johnson et al.l ( 20061 ) and is summarized below. 

The main part of our sample is comprised of 120 subgiant stars, which were 
selected from the Hipparcos catalog bas ed on the criteria 2.0 < My < 3.5, 0.55 < 
B -V < 1.0, and V < 7.6 (lESAlll997l fl Our sample of subgiants is illustrated 
in an H-R diagram shown in Figure [21 Also shown are the search domains of 
other Do ppler surveys containing evolved stars, along with the theoretical mass 
tracks of lcirardi et al.l (|2002), assuming solar abundances ([Fe/H]=0.0). 



Subgiants occupy an observational "sweet spot" in the H-R diagram. They 
exhibit relatively low levels of jitter, typically around 5 m s~^, which is only a 
factor of 2 higher than G dwarfs (Wright 2 005i) and signi ficantly lower than the 



20 m s ^ of jitter typical for giants (jHekker et al. 20061 ) . Like K giants, they 



have shed most of their primordial angular momentum and exhibit slow rotation 
velocities, with Vs'mi < 5 km s^^. Also, theoretical mass tracks along the 
subgiant branch are well separated, allowing for unambiguous mass estimates. 
Our sample of stars spans a stellar mass range 1.2 < M^/Mq < 2.2, which nearly 
doubles the stellar mass domain of the CCPS sample. 

Our planet search around subgiants has two primary goals. First, we wish 
to compare the orbital characteristics of planets around intermediate-mass stars 
to the large statistical ensemble of planets around lower-mass stars. Second, 
we wish to measure the fraction of stars with planets for stellar masses M* > 
1.3 Mq. To study the relationship between stellar mass and planet occurrence, 
we compare the planet fraction from our high-mass sample to that of the larger 
sample of FGK stars in the CCPS and the low-mass M dwarfs from the NASA 
Keck M Dwarf Planet Search (e.g. iButler et ani2006l ). 



^An additional 5 substellar companions have been detected around solar-mass giants 

^Our full survey contains an additional 39 giant stars with Mv < 2 and B — V < 0.85. However, 
these stars proved to be unsuitable Doppler targets, with a high fraction of close binaries and 
jitter in excess of 30 m s~^. We focus here only on the more stable subgiants. 
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Figure 2. H-R diagram illustrating the parameter space spanned by our 
sample of subgiants (circles). Also shown are the regions occupied by clump 
stars and K g iants (hashed regions); the theoretical stellar mass tracks of 
iGirardi et al.' (2002) for [Fe/H]=0.0 (thin lines); and the zero-age main se- 
quence (thick diagonal line). 



3. Planet Detections and Characteristics 



The fi rst detection from our sample of subgiants was announced in lJohnson et al 
(|2006i ): a short-period, Jovia n planet orbiting the 1.28 Mq subgiant HD 185269 
(see also Moutou et al. 20061 ) . With an orbital period P = 6.838 d and eccen- 



tricity e = 0.3, HD 185269 b has one of the highest eccentricities among the 
sample of known "hot Jupiters." The next batch of planets discovered from our 
sample orbit stars that are notably massive: HD 175541 (1.69 Mq), HD 192699 
(1.65 Mq) and HD 210702 (1.85 Mq, Johnson et al. 2007a). Following the the- 
oretical mass tracks of these three massive subgiants back to the main sequence 
reveals that they began life as early type dwarfs, with spectral types ranging 
from A5V to A8V. 

We have recently submitted for publication three additional long-period 
planet candidates orbiting the intermediate-mass subgiants n CrB (=HD 142091; 

= 1.80 Mq), HD 167042 (1.65 Mq), and HD 16175 (1.35 Mq; Johnson et al. 
2007, ApJ submitted). Figure [3] shows the distribution of semimajor axes and 
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Figure 3. Distribution of semimajor axes and eccentricities for planets 
around stars with M* < 1.5 Mq (filled circles) and subgiants with M* > 
1.5 M0 (open pentagrams). All of the planets around the evolved A-type 
stars orbit beyond ^ 1 AU. On the other hand, the eccentricity distributions 
for the two sets of planets are comparable. 



eccentricities of all known exoplanets orbiting stars with masses M* < 1.5 Mq. 
Also shown in the figure are planets orbiting subgiants with masses > 
1.5 Mq, including the 5 systems announced from our survey, our two strongest 
unpublished candi datej^, and two o ther planetary systems around massive sub- 
giants: HD 82744 (|Korzennik et al.ll2000l ) and HD 5319 (|Robinson et al.ll2007l ). 

While the eccentricities of planets around evolved A stars are very similar to 
those of planets around Sun-like stars. Figure [3] reveals a remarkable trend in the 
semimajor axes of planets around high-mass stars. Planets around evolved A 
stars (M» > 1.5 Mfr,) r eside preferentially in wide orbits at or beyond '^1.0 AU 
( Johnson et al. 2007bl ). This observed semimajor axis distribution of planets 



around high-mass stars differs significantly from that of planets around lower- 
mass stars, of which 51% orbit closer than 1 AU. This cannot be due to an 
observational bias, since Doppler shift measurements are most senstitive to giant 
planets in short-period orbits. While the radii of stars expand as they evolve 
away from the main sequence, the radii of subgiants are still small compared to 



■^Both unannounced systems have false alarm probabilities less than 1%, but lack sufficient phase 
coverage for publication at this time. 
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Figure 4. The fraction of stars with Jovian planets as a function of stellar 
mass. The error bars represent the uncertainties from Poisson statistics. The 
box above each bin shows the number of stars with detected planets, Nhosts^ 
and the total number of target stars, Nstars- 



the semimajor axis of even a P = 3 day hot Jupiter. Thus, it remains an open 
question as to whether the lack of close-in planets around A stars is related to 
the way planets formed in their nascent high-mass disks, or instead due to the 
effects of stellar mass on planet migration. 



4. Planet Occurrence vs. Stellar Mass 

Our sample of subgiants covers a range of stellar masses complementary to the 
CCPS sample of FGK stars and the sample o f low-m ass stars in the NASA Keck 
M Dwarfs Planet Search. In Johnson et al.l ( 2007al ) we showed that an analysis 



of the planet occurrence rate in three coarsely-spaced mass bins reveals a rising 
trend with stellar mass (Figure S]). For this analysis, we selected target stars 
and planet candidates with uniform detection characteristics, namely stars with 
more than 8 observations, and planets with a < 2.5 AU and Mp sini > 0.8 Mjup. 

The observed correlation between stellar mass and planet occurrence has 
important implications for planet formation theory, as well as for future planet 
search efforts. Stellar mass has now been identified as an additional sign post of 
planeticity, along with stellar metallicity, making A-type stars promising targets 
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for ground-based imaging surveys as well as space-borne astrometry and transit 
missions. 

5. Future Directions 

The study of planets orbiting massive stars is still in its infancy, with only 17 
systems currently known, compared to the 180 Sun-l ike and low-mass planet 
host stars discovered over the past decad^ ([Butler et al. 2006). Firmer conclu- 
sions about the occurrence rate and orbital properties of planets around A-type 
stars will require a much larger sample of detections. We have recently expanded 
our Doppler survey of subgiants to include 300 additional stars at Lick and Keck 
Observatories. Our primary goal is to confirm the correlation between stellar 
mass and planet occurrence seen in Figure [H If the ~ 9% occurrence rate holds, 
we expect to find 20-30 new planets over the next 3 years. This will represent 
a significant increase in the number of planets orbiting evolved A stars, and 
will allow us to perform a more robust analysis of the effects of stellar mass on 
planetary eccentricities, semimajor axes and multiplicity. 
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